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MODEL  STUDY  OF 

PORT  HURON  ICE  CONTROL  STRUCTURE 
Wind  stress  simulation 


D.S.  Sodhi,  D.J.  Calkins  and  D.S.  Deck 

INTRODUCTION 

One  function  of  an  ice  control  structure  (ICS) 
is  to  intercept  individual  ice  floes  in  the  early  winter 
and  retain  them  to  facilitate  formation  of  a  consol¬ 
idated  ice  cover.  The  ICS  to  be  installed  at  Port 
Huron,  Michigan,  will  serve  such  a  purpose,  forming 
a  consolidated  ice  cover  across  the  converging  shores 
of  Lake  Huron.  Already  these  shores  provide  support 
for  an  ice  arch  at  Port  Huron  and  hold  the  ice  back 
under  the  present  conditions  even  without  an  ICS. 

However,  the  naturally  formed  ice  arch  may  break 
due  to  the  anticipated  increase  in  winter  navigation, 
resulting  in  a  considerable  quantity  of  ice  passing 
into  the  St.  Clair  River  and  possibly  resulting  in  ice 
jams  in  its  lower  reach.  The  proposed  ICS  will  have 
an  opening  to  allow  ships  to  transit,  but  the  opening 
will  be  small  to  minimize  the  quantity  of  ice  going 
into  the  river. 

The  ice  cover  near  Port  Huron  consists  of  small 
floes  during  the  periods  of  freezeup  and  breakup  of 
the  ice  and  a  consolidated  ice  cover  for  the  rest  of 
the  winter.  In  midwinter,  winds  from  the  south 
sometimes  blow  the  ice  cover  to  the  north,  where 
the  ice  is  broken  up  by  the  wave  action;  reversal  of 
wind  direction  returns  the  ice  to  the  Port  Huron 
area  in  an  unconsolidated  state.  For  our  experiments 
an  unconsolidated  ice  cover  was  chosen  as  it  is  con¬ 
sidered  to  be  more  critical  in  terms  of  loading  the 
ICS  and  releasing  ice  through  the  opening.  Perham 
(1975)  reported,  from  force  data  taken  on  an  ice 
boom  in  the  Beauharnois  Canal,  that  force  levels 
were  high  during  the  freezeup  and  breakup  periods 
and  low  during  midwinter,  when  the  ice  cover  was 
fully  consolidated  from  shore  to  shore. 

The  model  study  described  in  this  report  was  de¬ 
signed  to  investigate  the  force  interaction  of  an  un¬ 
consolidated  ice  cover  with  the  shorelines  and  the 
ICS  at  Port  Huron,  and  also  to  determine  the  area 
of  ice  released  during  ship  transits.  An  experimental 
basin  was  built  in  which  the  uniform  wind  shear  stress 
on  an  ice  cover  was  simulated  by  water  flow  beneath 
the  ice.  The  forces  along  the  boundary  segments  of 
the  basin  and  the  ICS  were  measured  for  different 


directions  of  surface  shear  and  for  different  config¬ 
urations  of  the  ICS. 

This  report  presents  a  theoretical  model  for  de¬ 
termining  the  distribution  of  stresses  along  the  con¬ 
verging  boundaries  of  the  region.  This  theoretical 
distribution  is  compared  with  the  experimental  re¬ 
sults,  which  are  presented  in  dimensionless  form. 
Wind  data  (speed  and  direction)  at  Port  Huron  are 
summarized  for  the  years  1976-80,  and  estimates  are 
given  of  the  possible  forces  on  the  ICS  due  to  wind 
stresses  and  the  ice  release  through  the  opening. 


THEORETICAL  MODELS 

An  unconsolidated  ice  cover  consisting  of  a  large 
number  of  ice  floes  can  be  treated  as  a  continuum, 
provided  the  element  length  of  the  continuum  is 
large  compared  to  the  typical  floe  size  and  small  com¬ 
pared  to  the  scales  of  the  driving  forces  and  the  area 
under  investigation.  In  the  Arctic,  a  continuum  ele¬ 
ment  length  scale  on  the  order  of  100  km  has  been 
adopted  in  most  of  ice  dynamic  models  (Coon  1 980, 
Hibler  1980)  because  floe  sizes  rarely  exceed  10  km 
and  atmospheric  pressure  systems  are  typically  1000 
to  2000  km  across.  The  ratio  of  the  element  length 
scale  to  the  average  floe  size  is  taken  to  be  approxi¬ 
mately  1 : 10  in  the  modeling  of  ice  drift  in  the  Arctic. 

The  mass,  velocity  arid  forces  of  the  ice  for  these 
continuum  models  are  representative  averages  over  a 
region  large  enough  to  contain  many  ice  floes.  In  the 
following,  the  stress  in  the  ice  cover  is  represented 
by  the  stress  resultant  (N/m),  which  is  obtained  by 
integrating  the  stress  across  the  thickness  and  averaging 
the  stress  resultant  over  an  element  length.  The  use¬ 
fulness  of  a  continuum  approach  in  the  analysis  of 
the  stress  field  in  the  ice  cover  near  Port  Huron  will 
depend  on  the  size  of  ice  floes  relative  to  the  size  of 
the  ice  cover  area.  The  continuum  approach  is  justified 
during  freezeup  and  breakup  periods  as  the  region 
contains  a  large  number  of  floes.  Rigid-plastic  models, 
similar  to  the  models  for  granular  material,  have  been 
proposed  for  a  broken-up  ice  cover  (Pariset  and  Hauser 
1961,  Michel  1971,  Sodhi  1977,  and  Reimer  et  al. 
1979). 


Figure  7.  Geometry  of  an  ice  cover  held  between  converging  boun¬ 
daries  under  the  action  of  wind  stress  t.  The  coordinate  system  is 
represented  by  (x,  y)  or  (r,  0).  The  region  of  the  ice  cover  is  given 
r,  <  r  <  r2  and  -6W  <  $  <  dw.  The  components  of  the  stress  re¬ 
sultants  are  given  by  Nrr,  Nfl8  and  Nr9 . 


This  study  is  not  concerned  with  the  ice  transport 
in  the  lake.  Rumer  etal.  (1979)  have  presented  a 
general  summary  of  available  information  on  ice 
transport  in  the  Great  Lakes.  The  present  study 
deals  with  a  stationary  ice  cover  confined  against 
the  shorelines  by  the  wind,  though  small-scale  move¬ 
ment  and  deformation  in  the  form  of  ridging  and 
rafting  within  the  ice  cover  are  possible 

Figure  1  shows  the  geometry  of  an  unconsolidated 
ice  cover  being  contained  in  the  region  of  the  con¬ 
verging  shorelines  and  acted  upon  by  a  uniform  wind 
stress  t.  An  analogy  can  be  made  between  this  prob¬ 
lem  and  that  of  the  flow  and  arching  of  granular 
material  in  a  two-dimensional  wedge-shaped  hopper. 
The  analogy  exists  due  to  similarities  in  1 )  the  ice 
floes  and  grains  of  the  granular  material,  2)  the  sim¬ 
ilar  geometries  of  the  Port  Huron  shorelines  and  a 
converging  hopper,  and  3)  the  uniform  shear  in  the 
ice  cover  and  the  specific  weight  of  the  granular 
material.  Despite  of  these  similarities,  the  behavior 
of  an  ice  cover  and  that  of  a  confined  granular  ma¬ 
terial  have  several  differences  as  discussed  below. 

The  critical  condition  in  the  granular  material 
is  usually  described  by  a  Mohr-Coulomb  yield  cri- 


terion  which  is  characterized  by  two  parameters: 
cohesive  strength  (c)  and  the  angle  of  internal  fric¬ 
tion  (0).  In  addition,  it  is  implicitly  assumed  that 
the  deformation  of  the  material  is  in  shear  in  a  plane, 
i.e.  in  plane  strain  deformation.  Sodhi  (1977)  applied 
this  model  to  the  arching  and  drift  of  pack  ice  in  a 
wedge-shaped  channel  and  derived  an  expression  for 
the  critical  wind  stress  needed  to  cause  movement  of 
ice  in  the  channel.  Verification  of  the  usefulness  of 
the  model  was  attempted  by  making  comparisons 
with  ice  deformation  patterns  as  observed  with  satel¬ 
lite  imagery  in  the  Bering  Strait  region  and  in  the 
Amundsen  Gulf. 

In  a  similar  study,  Pritchard  et  al.  (1979)  use  a 
modified  yield  criterion  with  the  basic  characteristics 
of  a  Mohr-Coulomb  yield  criterion.  In  the  modified 
yield  criterion,  the  stress  resultants  are  always  nega¬ 
tive  and  of  more  than  a  certain  minimum  value  (-N*). 

These  modifications  in  the  yield  surface  are  justified 
because  an  ice  cover  composed  of  a  group  of  ice  : 

floes  cannot  support  any  tensile  stress  and  an  ice 
cover  can  resist  only  a  certain  maximum  level  of  com¬ 
pressive  stress  resultant  (-/V*),  beyond  which  it  de¬ 
forms  by  ridging,  rafting,  shore  pileup  or  rideup.  * 


Figure  2.  Yield  surface  of  an  unconsolidated  ice  cover. 


In  our  analysis,  the  yield  criterion  for  an  uncon¬ 
solidated  ice  cover  is  assumed  to  be  that  used  by 
Pritchard  et  al.  (1979).  It  is  depicted  by  a  Mohr 
circle  envelope  in  Figure  2a  and  by  a  yield  surface 
in  Figure  2b.  The  line  PQ  of  Figure  2a  may  be  de¬ 
fined  by  the  following  equation: 

(A/rr  +  Neg)  sin  0  -\(N„  -  Nee)*  +  4 /V^1'2 
=  2c  cos  0 

where  -N*  <  Nxr  <  0  and  -N*  <  Ngg  <  0.  If  the 
state  of  stress  point  in  the  domain  is  such  that  the 
Mohr  circle  is  tangent  to  line  PQ,  the  deformation  is 
in  the  shearing  mode  and  equivalent  to  the  principal 
stresses  being  represented  by  a  point  on  line  RS  or 
R'S'  of  Figure  2b.  When  the  Mohr  circle  is  tangent 
to  line  QQ',  the  deformation  is  in  the  compressional 
mode  and  equivalent  to  the  state  of  stress  on  lines 
ST  or  73’  of  Figure  2a.  The  state  of  stress  repre¬ 
sented  by  a  point  on  lines  OR  or  OR'  of  Figure  2b 
is  equivalent  to  the  Mohr  circle  being  tangent  to 
line  PPf  of  Figure  2a.  Point  R  in  Figure  2b  is  equiv¬ 


alent  to  the  unconfined  compressive  strength  (-/Va ) , 
which  is  an  important  parameter  related  to<he  co¬ 
hesive  strength  c  of  the  material,  as  given  in  Figure  2. 

In  the  plastic  analysis,  the  state  of  stress  is  re¬ 
stricted  to  lying  within  or  on  the  yield  surface.  To 
find  the  critical  state  of  an  ice  cover,  one  seeks  a 
case  in  which  the  stress  in  some  region  lies  on  the 
yield  surface,  and  in  order  to  simplify  the  analysis, 
we  will  consider  a  case  in  which  the  whole  region 
lies  on  the  vjeld  surface.  The  analysis  must  be  carried 
out  by  considering  the  different  cases  in  which  the 
state  of  stress  lies  on  the  three  individual  parts  of 
the  yield  curve,  i.e.  lines  OR,  RS  and  ST  of  Figure  2b. 

Case  1 

For  the  case  in  which  the  stress  resultants  in  the 
ice  cover  lie  on  line  RS  of  Figure  2b,  an  analysis  of 
the  stress  resultants  and  the  velocities  in  the  ice  cover 
is  given  in  Appendix  A.  The  distribution  of  stress 
resultants  along  the  boundary  of  the  converging  re¬ 
gion  is  also  derived  in  Appendix  A  for  a  particular  set 
of  boundary  conditions,  and  this  distribution  is  com¬ 
pared  with  the  experimental  results. 
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Case  2 

When  the  stress  resultants  lie  on  line  OR  of  the 
yield  surface  of  Figure  2b,  Reimer  et  al.  (1979)  have 
shown  that  only  materials  with  nonzero  unconfined 
compressive  strength  can  form  an  arch  to  withstand 
the  uniform  wind  'tress  and  other  external  loads  on 
the  arch. 

Case  3 

For  the  case  m  whu.h  the  ice  is  contained  ir,  a  re¬ 
gion  either  hs  jr  eontrol  structure  or  a  natural 
ice  arch,  u  is  possmle  tor  the  stresses  in  the  ice  cover 
to  be  in  such  a  slate  that  there  is  ridging,  rafting,  and 
shore  ice  piling.  In  such  a  situation,  the  state  of 
stress  in  the  ice  cover  is  represented  by  points  on  ST 
of  Figure  2b,  signifying  that  the  deformation  in  the 
ice  cover  would  tend  to  make  N2  (the  circumferen¬ 
tial  stress  component)  numerically  greater  than  A/, 

(the  radial  stress  component).  The  expressions  of 
N2  and  A/,  can  be  derived  in  terms  of  the  wind  stress 
t  and  the  limiting  compressive  stress  N*  as  proposed 
by  Reimer  et  al.  (1979).  In  a  favorable  situation, 
the  deformation  in  the  ice  cover  would  induce  an 
arch  to  form  across  the  converging  shorelines,  creating 
shore  pileups  and  a  rubble  field  relieving  the  ice 
pressure  on  the  ICS.  In  the  worst  situation,  it  is 
possible  that  the  ICS  would  be  loaded  by  the  limiting 
compressive  stress  resultant  N*  for  the  ice  cover. 

Though  several  theoretical  studies  have  been  con¬ 
ducted  to  determine  the  value  of  A'*  for  ridging, 
rafting  and  shore  pileups,  no  experimental  study  has 
been  conducted  to  verify  the  theoretical  results. 

Since  N*  depends  upon  such  processes  as  ridging, 
rafting,  ice  pileups  on  shores  and  instability  of  ice 
blocks,  it  is  difficult  to  determine  its  value  for  a 
given  ice  cover,  but  Prichard  etal.  (1979)  have  used 
a  value  of  1.5  x  105  N  m'1  for  N*  in  their  analysis 
of  sea  ice  flow  through  straits. 

EXPERIMENTAL  PROGRAM 
Experimental  facility 

The  simulation  of  wind  stress  on  the  ice  cover 
over  a  large  hydraulic  (Froude)  model  area  is  diffi¬ 
cult,  because  of  the  problems  of  blowing  air  uni¬ 
formly  over  a  large  area  and  of  maintaining  a  velocity 
gradient  in  the  boundary  layer  to  impart  the  required 
magnitude  of  shear  stress  on  the  ice  cover.  Since 
these  difficulties  could  not  be  overcome,  it  was  pro¬ 
posed  that  a  wind  stress  simulation  model  should 
be  constructed  separate. y  from  the  hydraulic  model. 

In  this  model,  a  uniform  shear  stress  on  the  floating 
ice  cover  is  imparted  by  flowing  water  beneath  the 
ice  instead  of  blowing  air  over  it. 


A  large  basin,  approximately  15.2  x  15.9  m  (50 
x  52  ft),  was  constructed  of  wood  and  sealed  with 
fiberglass  resin.  The  water  level  was  controlled  with 
a  gate  at  the  downstream  end  of  the  model.  The 
basin  could  handle  0.34  m3/s  (12  ft3/s),  of  water, 
which  resulted  in  a  water  velocity  of  0.15  m/s  (0.5 
ft/s)  at  a  depth  of  about  0.14  m  (0.46  ft).  The  veloc¬ 
ity  of  the  water  could  be  increased  by  maintaining 
a  lower  water  depth  in  the  basin,  but  this  resulted 
in  instability  and  underturning  of  ice  floes.  Usually 
a  flow  rate  of  0.32  m3/s  (1 1 .3  ft3/s)  was  maintained 
to  have  a  depth  of  0.18  m  (0.58  ft)  and  an  average 
water  velocity  of  0.1 1  m/s  (0.37  ft/s)  under  the  ice 
cover. 

Figure  3  shows  the  general  layout  of  the  basin 
area  and  the  boundary  elements  which  skimme'4  the 
water  sui  face  from  above  and  held  back  the  rru  Jel 
ice  from  drifting  downstream.  The  0.5-m-long 
boundary  elements  were  attached  to  instrumented 
rods  hung  from  a  square  aluminum  support  beam. 

A  photograph  of  the  boundary  elements,  instru¬ 
mented  rods  and  support  beams  is  shown  in  Figure  4. 
The  6-m-  (19.7-ft-)  long  unit,  containing  12  boundary 
elements,  could  be  placed  at  different  locations  rel¬ 
ative  to  the  direction  of  water  flow  to  simulate  wind 
from  different  directions. 

The  instrumented  rods,  capable  of  measuring 
forces  in  two  directions,  were  connected  to  a  data 
acquisition  system  controlled  by  a  computer.  The 
data  acquisition  system  was  housed  in  a  temperature- 
controlled  box.  The  instrumented  rods  were  cali¬ 
brated  individually  and  the  calibration  matrices  of 
all  rods  were  stored  on  a  tape.  At  the  time  of  taking 
the  force  measurements,  98  channels  were  scanned 
simultaneously,  with  the  readings  transformed  to 
force  data  and  the  averages  of  eight  readings  recorded. 

Although  the  force  measurement  system  had  a 
resolution  of  0.01  N,  the  system  was  sensitive  to 
temperature  variations.  Considerable  care  had  to  be 
taken  to  ensure  that  there  was  little  or  no  drift  in 
the  force  measurements.  Initial  and  final  readings  at 
the  beginning  and  end  of  the  experiments  were  com¬ 
pared  to  determine  if  zero  drift  had  occurred  in  the 
force  measurement. 

Three  types  of  model  ice  were  used:  1 )  polyethy¬ 
lene  plastic  pieces  10  x  10  x  0.6  cm  (4  x  4  x  1/4  in.), 

2)  freshwater  ice  (randomly  shaped  ice  floes),  and 

3)  urea  ice  (randomly  shaped  ice  floes).  The  fresh¬ 
water  or  urea  ice  sheet  was  grown  in  still  water,  and 
it  was  broken  up  to  form  a  field  of  floes  of  random 
shape  and  given  average  size.  After  breaking  the  ice 
sheet,  the  water  flow  in  the  basin  was  started,  and 
the  water  level  was  adjusted  to  the  level  of  the  boun¬ 
dary  elements  which  skimmed  the  water  surface  and 
kept  the  ice  in  the  region  upstream  of  them.  The  ice 
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Figure  3.  Sketch  of  the  test  basin. 


downstream  from  the  boundary  elements  was  allowed 
to  flow  over  the  gate  into  the  storage  tank  where  it 
melted  and  kept  the  water  temperature  close  to  0°C. 

The  shear  stress  imparted  by  the  water  to  the 
model  ice  cover  was  determined  by  summing  the 
measured  force  components  in  the  streamwise  direc¬ 
tion.  The  total  force  was  then  divided  by  the  area 
of  the  cover  to  obtain  the  uniform  shear  stress  t 
acting  on  the  ice  cover. 

Scaling  factors 

In  order  to  relate  model  results  to  prototype 
values,  the  experimental  results  are  presented  in  a 
nondimensional  form  that  can  be  used  for  both  the 
model  and  the  prototype.  The  expressions  derived 
for  the  stress  resultants  in  Appendix  A  have  been 
normalized  with  respect  to  the  product  of  shear  stress 
r  and  the  distance  d  of  the  ice  cover  region  shown 
in  Figure  5.  For  the  corresponding  points  in  the 
prototype  and  the  model,  the  normalized  stress  com¬ 
ponents  are  functions  of  the  position  of  the  point, 
provided  the  behavior  of  the  broken-up  model  ice 
cover  represents  the  prototype  behavior.  Referring 
to  Figure  5,  we  have 


=  f  (position  of  a  point). 


Figure  5.  Scale  factors  for  the  wind  stress 
simulation. 


One  can  use  the  above  equation  for  the  prototype 
condition  in  the  following  manner: 

N  ~  (N*\ 

^ij^prototype  ^T<^prototype  \Td) 

'  '  model 

The  distribution  of  normalized  stress  (N^/rd) 
is  determined  from  model  experiments.  The  factor 
(rd)  for  the  prototype  is  estimated  from  the  wind 
data  and  the  extent  of  the  ice  cover. 

Experimental  results 

The  results  of  the  experiments  are  presented  in 
three  categories:  1 )  distribution  of  stress  components 
along  the  converging  boundaries  without  an  ICS,  2) 
normalized  loads  for  various  configurations  of  the 
ICS  and  3)  the  area  of  ice  released  during  ship  tran¬ 
sits. 

Distribution  of  stress  components  along  the  con¬ 
verging  boundaries 

The  region  of  the  floating  ice  cover  for  these  ex¬ 
periments  is  shown  in  Figure  1.  The  included  angle 
between  the  boundaries  is  84°,  approximately  the 
angle  between  the  shore-fast  ice  edges  at  Port  Fluron. 
The  direction  of  the  uniform  shear  stress  is  parallel 
to  the  line  of  symmetry  of  the  region.  Both  plastic 
pieces  and  freshwater  ice  were  used  in  the  experiments 
to  determine  any  significant  differences  in  their  be¬ 
havior. 

The  main  reason  for  conducting  this  set  of  ex¬ 
periments  was  to  verify  the  theoretical  distribution 
of  normal  stress  resultants  along  the  boundary  as 
given  by  eq  A26.  The  experiments  were  conducted 
so  that  the  ice  cover  would  move  upstream  a  small 
distance  from  the  boundaries.  Then  it  would  travel 
forward  and  come  to  rest  against  the  boundary  ele¬ 
ments  after  deformation  took  place  within  the  ice 
cover,  due  to  the  interaction  of  ice  with  the  boun¬ 
daries.  The  deformation  was  equivalent  to  creating 
a  critical  state  of  stress  in  the  ice  cover,  and  the  dis¬ 
tribution  of  stress  along  the  boundary  was  fixed  once 
the  ice  cover  stopped  moving. 

Figure  6a  is  a  plot,  divided  into  0.5-m-long  boun¬ 
dary  segments,  of  the  typical  force  measurement 
across  the  basin  in  the  normal  and  tangential  direc¬ 
tions  relative  to  the  boundary.  The  distribution  of 
the  stress  components  shows  a  trend,  but  there  is 
considerable  variation  from  segment  to  segment. 

This  variation  results  from  the  measuring  elements 
not  being  long  enough  compared  to  the  floe  size  to 
give  an  average  representative  value  of  the  stress  re¬ 
sultant.  When  the  same  data  are  plotted  in  Figure  6b 
by  summing  stress  resultants  over  1-m-long  boundary 
segments,  the  trend  of  representative  stress  resultants 
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a.  Measuring  element  length  =  0.5  m,  freshwater  ice  floes  used. 


Figure  6.  Distribution  of  measured  boundary  stress  in  the  normal 
and  tangential  directions  and  the  theoretical  distribution  of  normal 
stresses  (eg  A  26,  <j>  -  20°), 
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Figure  9.  Same  data  presentation  as  in  Figure  6b  for  polyethy¬ 
lene  pieces  (10  x  70  x  0.6  cm). 


can  be  easily  seen.  Figure  7  shows  the  same  type 
of  plot  as  Figure  6b  but  for  a  force  measurement 
taken  4  minutes  later.  The  ratio  of  the  measuring 
length  of  (1  m)  to  the  average  floe  size  (10  cm)  is 
approximately  10,  a  ratio  used  in  many  problems 
that  are  concerned  with  flow  in  granular  media. 

The  theoretical  distribution  of  stress  resultants 
at  the  boundary  is  given  in  eq  A26  for  a  friction  less 
boundary  condition.  The  measured  tangential  forces 
along  the  boundary  are  low,  and  thus  comparison 
of  the  theoretical  and  experimental  distributions 
of  the  normal  stress  resultant  at  the  boundary  is 
justified.  Figures  6  and  7  also  show  the  distribution 
of  normal  stress  components  at  the  boundary  as 
given  by  eq  A26.  Similar  distributions  of  theoretical 
and  measured  normal  stress  components  arc  pre¬ 
sented  in  Figures  8  and  9  for  freshwater  and  poly¬ 
ethylene  pieces,  respectively. 

The  similarity  between  the  theoretical  and  meas¬ 
ured  distribution  of  the  normal  stress  resultant  is 
good  except  at  a  few  points.  From  these  compari¬ 
sons,  it  is  deduced  that  a  fragmented  ice  cover  can 
be  modeled  with  a  Mohr-Coulomb  yield  surface  and 
that  the  angle  of  internal  friction  is  about  20°  for 
freshwater  ice  and  is  about  10°  for  plastic  pieces. 

Forces  on  the  ice  control  structure 

Experiments  were  conducted  to  determine  the 
forces  exerted  by  the  ice  cover  on  the  ICS,  which 
was  attached  to  transducer  rods  at  the  downstream 


end  of  the  ice  cover  as  shown  in  Figure  10.  The 
measured  forces  on  each  end  of  the  ICS  were  added 
to  determine  the  force  per  unit  length  of  the  ICS, 
which  was  then  normalized  by  rd  (see  Fig.  10). 

For  each  test,  a  particular  configuration  of  the 
ICS  was  chosen,  and  the  force  measurements  were 
taken  after  disturbing  the  ice  cover.  This  procedure 
was  repeated  several  times  during  the  course  of  a 
test,  and  the  average  normalized  forces  per  unit 
length  of  the  ICS  were  computed  for  each  test.  These 
average  values  for  all  the  tests  are  presented  in  Tables 
1  and  2  for  plastic  pieces  and  freshwater  ice,  re¬ 
spectively. 

Many  problems  were  encountered  during  measure¬ 
ment  of  forces  due  to  humidity  and  temperature 
variations  in  the  instrument  box.  The  results  of  the 
ice  force  tests  do  not  indicate  clearly  any  particular 
pattern,  but  the  following  two  conclusions  may  be 
made  with  caution: 

1 .  As  the  distance  from  the  apex  to  the  ICS 
(distance  e  in  Fig.  10)  decreases,  the  ICS  seems 
to  move  progressively  under  the  protection  of 
an  arch,  and  thus  experiences  smaller  forces 
from  the  ice  cover. 

2.  The  orientation  of  the  ICS  affects  the  force 
transmitted  to  it  by  assisting  in  the  formation 
of  an  ice  arch  across  the  converging  boundaries. 

As  mentioned  earlier,  these  conclusions  are  not 
clearly  indicated  by  the  results  of  the  small  number 
of  tests,  but  there  is  a  trend  towards  them  in  the 
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Figure  10.  Geometry  of  an  ice  control  structure  for  force  and 
ice  release  measurements. 


Table  1.  Average  normalized  total  force 
per  unit  length  of  the  ice  control  struc¬ 
ture  when  using  polyethylene  pieces 
(10  x  10  x  0.6  cm). 

A  verage 


Experiment 

a,  0 

e 

(m) 

SI 

(m) 

value  of 
N/rd 

ii 

0° 

1.10 

0.60 

0.98 

12 

23° 

1.10 

0.70 

0.95 

13 

51° 

1.10 

1.00 

0.69 

15 

0° 

2.00 

1.20 

0.85 

16 

0° 

1.10 

0.60 

1.23 

18 

0° 

0.76 

0.30 

0.71 

results.  More  comments  will  be  made  on  these  re¬ 
sults  in  the  Summary  and  Conclusions. 

Area  of  ice  released  during  ship  transits 

Experiments  were  conducted  to  measure  the 
area  of  ice  released  through  the  opening  in  the  ICS 
when  the  ice  cover  was  disturbed  by  moving  a  blunt 
object,  representing  ship  transits  through  the  ice 
cover.  The  tests  were  conducted  for  different  orien¬ 
tations  (a,  0)  and  positions  of  the  ICS  (e)  and  for 
different  sizes  of  the  opening  (see  Fig.  10).  These 
tests  were  conducted  using  plastic  pieces  and  broken- 
up  freshwater  and  urea  ice.  The  average  flow  size 


Table  2.  Average  normalized  total  force 
per  unit  length  of  the  ice  control  struc¬ 
ture  when  using  freshwater  ice. 


tiperiment 

a,  |  j 

c 

(m) 

C 

(m) 

A  i ’cragc 

value  of 
N/rd 

60 

34° 

3.65 

2.74 

1.16 

61 

34“ 

3.65 

2.74 

1.15 

62 

0° 

3.65 

1.83 

2.22 

63 

-30° 

3.65 

2.59 

0.31 

64 

-30° 

3.65 

2.59 

0.67 

65 

36° 

3.65 

2.44 

- 

66 

36° 

3.65 

2.44 

2.57 

67 

36° 

3.65 

2.44 

2.52 

68 

36° 

3.65 

2.44 

1.76 

69 

45° 

1.98 

1.37 

1.46 

70 

0° 

1.98 

1.22 

0.63 

71 

0° 

1.98 

1.22 

1.94 

72 

0° 

1.98 

1.22 

0.82 

73 

0° 

1.98 

1.22 

0.62 

74 

-36° 

1.98 

1.52 

0.71 

IS 

45° 

1.98 

1.37 

1.29 

76 

45° 

1.98 

1.37 

- 

77 

70° 

1.98 

1.52 

1.26 

a  was  determined  by  weighing  a  number  of  floes. 
The  area  of  ice  released  during  each  ship  transit  was 
determined  by  weighing  the  ice.  The  area  of  ice  was 
then  normalized  with  respect  to  the  square  of  the 
opening  in  the  ICS,  and  the  mean  value  of  the  nor¬ 
malized  area  of  ice  released  was  determined  for  20 
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Figure  7  7.  Plot  of  mean  normalized  ice  discharge  A,/b2  with  respect 
to  a/b  for  upbound  ship  passages. 


Figure  12.  Plot  Of  mean  normalized  ice  discharge  Ar/h2 
with  respect  to  a/b  for  downbound  ship  passages. 

upbound  and  20  downbound  passages  of  ship.  The 
ice  release  during  the  ship  transits  was  a  random  phe¬ 
nomenon  that  depended  upon  the  arrangement  of 
ice  floes  near  the  opening  before  the  occurrence  of 
arching.  The  ability  to  arch  has  been  demonstrated 
to  be  a  function  of  b/a  by  Calkins  and  Ashton  (1975). 


The  mean  normalized  ice  area  released  ( AJb 2) 
during  ice  release  experiments  is  plotted  with  respect 
to  the  ratio  of  the  average  floe  size  to  the  opening  in 
the  ICS  (a/b)  in  Figures  1 1  and  1 2  for  the  upbound 
and  downbound  ship  passages,  respectively.  The 
symbols  in  Figures  1 1  and  12  refer  to  the  results  of 
a  particular  test  conducted  for  a  particular  configur¬ 
ation  of  the  ICS.  The  details  of  each  particular  test 
are  given  in  Table  3. 

No  definite  trend  can  be  seen  in  the  results  of 
these  tests  since  the  data  are  so  scattered.  The  values 
of  mean  normalized  ice  area  released  are  higher  than 
those  obtained  from  the  tests  conducted  in  the  hy¬ 
draulic  model  (Calkins  ct  al.  in  prep).  This  increase 
is  perhaps  attribu  able  to  high  values  of  shear  stress 
t  acting  on  the  ice  cover,  as  they  are  an  order  of  mag¬ 
nitude  higher  than  those  for  the  tests  conducted  in 
the  hydraulic  model.  The  only  conclusions  that  can 
be  reached  are  that  1 )  the  ice  cover  arches  for  low 
values  of  a/b  (0.1 1),  and  2)  the  mean  values  of  nor¬ 
malized  ice  area  released  (Ar/b2)  are  less  than  3  for 
downbound  ship  passages  and  less  than  2  for  upbound 
ship  passages,  with  the  exception  of  the  data  from 
one  particular  test. 

The  results  of  the  ice  release  data  (Table  3)  indi¬ 
cate  less  release  of  ice  when  the  wind  is  blowing  from 
a  direction  other  than  along  the  centerline.  The  de¬ 
crease  in  ice  release  is  associated  with  the  oblique  in¬ 
cidence  of  the  opening  in  the  ICS  relative  to  the  ice 
floe  movement. 


Table  3.  Results  of  ice  release  experiments*. 


t  x  peri  mem 
number 

Type  of 
icef 

Orientation 
of  ICS 

Distance 
from  apex 

t 

(m) 

Floe  size 

a 

(mm) 

Size  of 
opening 
b 

(cm) 

fee 

thickness 

t 

(mm) 

Symbol  ft 

Ratio 
of  floe 
to  opening 

ilb 

Ice  release 

Ar/b2 

Upbound  Downbound 

Wind 

Direction 

0 

60 

W 

34° 

3.65 

123 

51 

a 

0.24 

0.54 

. 

0° 

61 

w 

34° 

3.65 

159 

51 

_ 

b 

0.31 

0.43 

- 

of 

62 

w 

0° 

3.65 

115 

51 

11 

0.23 

0.98 

- 

0° 

63 

w 

-30“ 

3.65 

119 

51 

15 

d 

0.23 

1.04 

- 

0° 

64 

w 

-38“ 

3.65 

92 

5! 

30 

e 

0.18 

1.03 

- 

0° 

65 

w 

36° 

3.65 

86 

51 

18 

f 

0.17 

0.76 

- 

0° 

66 

w 

36° 

3.65 

80 

51 

11 

g 

0.16 

0.39 

- 

0° 

67 

w 

36“ 

3.65 

91 

51 

9 

h 

0.18 

1.56 

- 

0“ 

68 

w 

36“ 

3.65 

113 

51 

25 

i 

0.22 

1.29 

1.46 

0° 

69 

w 

45° 

1.98 

88 

51 

13 

0.17 

0.40 

1.63 

0“ 

70 

w 

0“ 

1.98 

97 

51 

23 

k 

0.19 

1.67 

2.31 

0° 

71 

w 

0“ 

1.98 

11 6 

51 

26 

1 

0.23 

1.08 

2.05 

0” 

72 

w 

0° 

1.98 

122 

51 

14 

m 

0.24 

0.47 

1.96 

0“ 

73 

w 

0° 

1.98 

110 

51 

13 

n 

0.22 

0.68 

2.17 

0° 

74 

w 

-36“ 

1.98 

83 

51 

15 

0 

0.16 

0.6 

2.3 

0° 

75 

w 

45° 

1.98 

79 

34 

13 

P 

0.23 

0.47 

2.28 

0“ 

76 

w 

45“ 

1.98 

75 

24 

9 

q 

0.31 

0.78 

2.04 

0° 

77 

w 

45“ 

1.98 

87 

80 

13 

r 

0.11 

0.14 

0.72 

0° 

78 

u 

45° 

1.98 

96 

80 

12 

s 

0.12 

0.31 

1.05 

0“ 

79 

u 

45° 

1.98 

83 

51 

21 

t 

0.16 

3.14 

3.51 

0“ 

80 

w 

45” 

1.98 

95 

51 

15 

u 

0.19 

0.58 

0.97 

0° 

81 

w 

45° 

** 

113 

41 

15 

V 

0.28 

1.38 

2.70 

0° 

82 

w 

45° 

*• 

82 

41 

26 

w 

0.20 

0.69 

1.63 

0° 

83 

w 

45° 

•  « 

117 

51 

18 

X 

0.16 

0.49 

1.55 

0° 

90 

w 

45° 

1.98 

100 

61 

20 

y 

0.16 

1.42 

1.61 

21° 

91 

w 

45° 

1.98 

94 

41 

19 

2 

0.23 

0.08 

2.33 

21° 

92 

w 

45° 

1.98 

126 

51 

23 

A 

0.25 

0.69 

2.26 

21° 

93 

w 

45“ 

1.98 

89 

51 

25 

B 

0.17 

0.05 

1.17 

42° 

94 

w 

45° 

1.98 

95 

61 

25 

C 

0.16 

0.50 

1.28 

42° 

95 

w 

45° 

1.98 

93 

81 

24 

D 

0.11 

0.23 

0.86 

42“ 

96 

w 

45“ 

1.98 

217 

81 

25 

E 

0.27 

0.20 

1.30 

42“ 

27 

p 

45° 

4.00 

100 

48 

6 

F 

0.21 

1.13 

_ 

0“ 

29 

p 

45” 

4.00 

100 

33 

6 

G 

0.31 

0.84 

. 

0° 

30 

p 

45” 

4.00 

100 

25 

6 

H 

0.40 

0.48 

- 

0° 

31 

p 

0” 

4.00 

100 

25 

6 

1 

0.40 

0.6 

- 

0° 

32 

p 

0” 

4.00 

100 

35 

6 

i 

0.31 

0.86 

- 

0° 

33 

p 

0° 

4.00 

100 

48 

6 

K 

0.21 

1.02 

- 

0° 

34 

p 

-30° 

4.00 

100 

48 

6 

L 

0.40 

0.84 

- 

0“ 

40 

p 

-30° 

4.00 

100 

33 

6 

M 

0.21 

0.64 

- 

0“ 

41 

p 

•30° 

4.00 

100 

25 

6 

N 

0.31 

0.64 

- 

0* 

43 

p 

45° 

4.00 

100 

48 

6 

O 

0.21 

1.24 

_ 

0° 

45 

p ♦** 

45° 

4.00 

200 

102 

6 

P 

0.20 

0.72 

- 

0“ 

113 

p 

45“ 

4.00 

100 

51 

6 

0 

0.20 

1.26 

- 

30° 

114 

p 

45° 

4.00 

100 

51 

6 

R 

0.20 

0.80 

“ 

45° 

*  Refer  to  F igure  10  for  the  interpretation  of  each  letter  listed  in  this  table, 
t  W-freshwater  ice,  P-plastic  pieces,  U-urea  ice. 

**  Double  booms  were  used  for  these  tests  as  shown  in  Figure  10. 
ft  The  following  symbols  are  used  to  plot  the  results  in  Figures  1 1  and  12. 
***  Randomly  shaped  plastic  pieces. 


ANALYSIS  OF  WIND  DATA 
FOR  LOWER  LAKE  HURON 

The  data  on  wind  conditions  in  the  study  area 
were  obtained  from  three  sources.  The  U.S.  Army 
Detroit  District,  Corps  of  Engineers  received  a  wind 
analysis  for  January,  February  and  March,  for  Sarnia, 
Ontario  (from  Acres  American  Inc.  1974).  The 
National  Oceanic  and  Atmospheric  Administration 
(NOAA)  published  a  regional  set  of  meteorological 
parameters  for  Southern  Lake  Huron  for  1960-73, 
but  this  summary  omits  January  through  March. 


The  third  data  set  comes  from  the  U.S.  Coast  Guard 
Station  at  Port  Huron,  and  is  on  file  for  the  years 
1975-1980. 

The  Acres  American  Inc.  (1974)  study  showed 
that  the  mean  maximum  hourly  wind  speed  was 
roughly  20  mph  from  all  eight  quadrants.  The  max¬ 
imum  wind  speed  was  39  mph  from  the  SW  quad¬ 
rant,  and  this  was  defined  as  the  maximum  one-hour 
value.  The  NOAA  data  indicated  that,  for  the  0.1% 
event  between  1960-73,  the  wind  speed  was  greater 
than  48  mph,  but  only  for  those  nine  months  indi¬ 
cated  in  Table  4. 


Table  4  .  Summary  of  wind  data  (Acres  American,  Inc.  1974  and 
NOAA  1960-73). 


Mean  wind  speed  (mph)  Maximum 

jan-Mar  Apr-Dee  Frequency  Frequency  wind  %peed 


Wind  (1960-66)  (1960-73)  (%)  (%)  ) an-Mar 

direction  Acres _ NOAA _ Apr-Dee  /  an-Mar _ (mph) 


N 

11.7 

16.4 

12.2 

5 

32 

NE 

10.4 

14.7 

9.8 

9 

25 

E 

9.0 

13.4 

8.1 

7 

29 

SE 

10.8 

14.0 

11.6 

IS 

33 

S 

10.1 

15.1 

16.0 

11 

30 

SW 

11.6 

15.9 

15.2 

16 

39 

W 

10.8 

16.9 

14.6 

19 

36 

NW 

11.6 

16.1 

11.8 

18 

33 

Table  5.  Fort  Gratiot  maximum  recorded  wind  velocities 

(mph)  (maximum  3-hr  velocity/direction). 

1974  -  75  1  975-76 

1976-77 

1977-78  1  978-79 

1979-80 

Dec 

* 

32/SW 

37/NE 

30/SW 

30/SW 

|an 

35/S  27/N 

28/NE 

30/NE 

27/N 

35/SE 

Feb 

30/SW  30/E 

25/NW 

27/N 

20/NW* 

33/N 

Mar 

30/E  25/SW 

40/NW 

26/N 

30/SE 

40/ SW 

♦Partial  record  for  that  month. 


The  frequency  of  duration  is  obviously  different 
depending  upon  the  months  considered.  The  data 
in  Table  4  show  that  the  wind  direction  has  a  pre¬ 
ferred  SW  to  NW  pattern  during  the  winter  months. 
The  wind  direction  is  uniformly  distributed  during 
the  other  nine  months.  The  mean  wind  speed  is 
slightly  higher  for  the  nonwinter  months. 

Data  from  the  Fort  Gratiot  Coast  Guard  Station 
at  Port  Huron  have  been  reported  for  6-hr  quarters 
with  the  basic  data  collected  at  3-hr  intervals.  The 
Detroit  District  abstracted  the  information  from 
their  files  to  a  computer  base  file.  The  data  given  in 
Table  5  are  abstracted  from  the  basic  data  and  repre¬ 
sent  the  maximum  monthly  3-hr  wind  velocity. 

Table  5  summarizes  the  frequency  of  wind  duration 
for  the  five  years  on  the  monthly  basis.  The  months 
of  January  and  February  show  a  strong  SW  to  NW 
pattern  of  wind  direction  while  March  has  a  more 
balanced  pattern.  The  monthly  summaries  for  five 
years  (1976-80)  are  given  in  Appendix  B.  A  maximum 
wind  speed  of  30  to  40  mph  occurs  occasionally. 

Winds  of  shorter  duration,  but  higher  speed,  have 
probably  not  been  included  in  these  records.  The 
six  years  of  wind  data  from  Fort  Gratiot  show  the 
same  general  trend  as  that  found  by  Acres  American, 
Inc.  (1974).  The  wind  directions  from  the  southwest 
to  northwest  quadrants  account  for  more  than  50% 
of  the  duration. 


SUMMARY  AND  CONCLUSIONS 

The  first  set  of  experiments  was  conducted  to 
determine  the  distribution  of  stress  components  at 
the  boundaries  of  a  converging  region  of  ice  cover 
acted  upon  by  wind  stress,  and  the  results  are  pre¬ 
sented  in  nondimensional  form.  Hence  these  results 
can  be  applied  to  a  region  of  any  size  for  which  the 
behavior  of  an  ice  cover  can  be  represented  by  that 
in  the  model  experiments. 

In  the  second  set  of  experiments,  an  ice  control 
structure  was  installed  in  the  model  with  the  objec¬ 
tive  of  determining  the  amount  of  ice  released  through 
the  opening  in  the  ICS  and  the  ice  forces  on  the  ICS. 
In  this  case,  a  length  scaling  factor  can  be  obtained 
from  the  similarity  of  the  prototype  and  model 
geometry. 

Release  of  ice  through  the  opening  of  an  ICS 

Referring  to  Table  3  or  Figures  1 1  and  12,  wc 
conclude  that  the  ice  cover  arches  for  low  values  of 
alb  a  0.1 1  to  0.1 5,  and  that  the  mean  average  value 
of  normalized  ice  area  (AT/b2)  released  is  less  than  3. 

If  the  average  size  of  prototype  ice  floes  is  about 
20  m  (65  ft),  the  ice  cover  would  arch  across  an 
1 33-m-  (437-ft-)  wide  opening  in  the  ICS  for  a/b  = 

0. 1 5.  If  the  size  of  the  opening  is  selected  to  76  m 
(250  ft)  and  a/b  =  0.1 5,  an  ice  cover  will  arch  when 
the  average  size  of  the  floes  is  1 1 .4  m  (37  ft). 
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Taking  a  figure  of  AJb-  =  3,  an  estimate  of  icc 
release  for  100  round  trips  of  ships  through  an  open¬ 
ing  of  1 22  m  (400  ft)  in  the  ICS  is  approximately 
9  km2  (3.45  mile2).  This  figure  should  be  multi¬ 
plied  by  a  factor  for  the  probability  of  ship  transits 
during  high  winds. 

Ice  forces  on  the  ice  control  structure 

While  conducting  ice  release  experiments,  we 
measured  the  ice  forces  on  the  model  ICS,  and  the 
normalised  icc  forces  per  unit  length  of  ICS  arc 
given  in  Table  1  and  2.  From  Figure  10  and  the 
probable  location  of  the  ICS  near  Port  Huron,  a 
length-scaling  factor  can  be  determined  by  taking 
the  ratio  of  distances  («?)  in  the  model  and  in  the 
prototype.  The  approximate  location  of  the  ICS 
along  the  U.S.  and  the  Canadian  shorelines  is  known 
to  be  762  m  (2500  ft)  uplake  from  the  mouth  of 
the  St.  Clair  River.  The  corresponding  lengths  in 
the  model  are  1 .98  m  (6.50  ft)  and  3.65  m  ( 1 2  ft) 

(sec  Table  2  and  Fjg.  10).  Thus  the  length  scale 
factor  X  is  equal  to  385  or  209,  depending  upon  the 
distances  (e)  in  the  model.  For  the  case  of  X  =  385, 
a  cover  length  of  1 1 .50  m  (37.75  ft)  in  the  model 
corresponds  to  4.43  km  (2.75  miles)  of  cover  length 
in  the  prototype.  This  cover  length  is  not  long  enough 
to  be  certain  that  an  icc  cover  beyond  this  distance 
would  not  load  the  ICS  more  than  the  ice  cover  with¬ 
in  this  distance.  A  longer  icc  cover  length  could  be 
achieved  cither  by  having  a  larger  model  area  or  by 
conducting  a  test  with  larger  values  of  the  length 
scale  factor  X.  The  first  approach  would  have  meant 
constructing  a  larger  test  facility.  The  second  approach 
has  limitations  in  terms  of  conducting  a  test  with  a 
very  small  model  of  the  ICS  and  even  smaller  ice 
floes.  For  the  ice  release  experiments,  the  ice  cover 
was  sufficiently  large  because  icc  release  occurs  from 
an  area  adjacent  to  the  ICS. 

Despite  the  insufficient  length  of  the  model  ice 
cover,  it  is  useful  to  estimate  the  ice  forces  per  unit 
length  of  the  ICS  from  the  force  measurements 
taken  during  the  ice  bleeding  experiments.  The 
values  of  normalized  total  load  (N/rd)  on  the  ICS 
(given  in  Table  3)  range  from  0.31  to  2.52,  and  their 
average  value  is  1 .3.  The  variation  in  these  measure¬ 
ments  is  attributed  to  the  possible  problems  of  drift 
in  the  force-measuring  system  due  to  temperature 
variations,  although  great  care  was  taken  to  avoid 
them.  The  variations  may  also  result  from  different 
arrangements  of  ice  floes  and  changes  in  the  behavior 
of  the  icc  cover.  Despite  the  variation,  the  normal¬ 
ized  force  (N/rd)  on  the  ICS  is  taken  to  be  equal  to 
1 .0  in  the  following  discussion. 

From  the  wind  data  for  Port  Huron,  a  maximum 
wind  velocity  of  about  1 8  m/s  (40  mph)  is  anticipated. 


The  value  of  wind  shear  stress  is  estimated  from  the 
expression,  r  =  cdpV2 ,  in  which  a  value  of  0.003 
for  cd  and  1.3  kg/m3  for  p  (density  of  air)  is  taken. 
Thus,  t  =  1 .26  N/m2  (1 ,83x  10‘4  psi)  for  a  wind 
velocity  of  18  m/s  (»40  mph).  For  a  cover  length 
of  4.43  km  (2.75  miles),  we  get  a  value  of  (rd) pr0|0,vpe 
=  5600  N/M  (383  Ibf/ft).  If  we  assume  the  thickness 
of  ice  to  be  0.5  m,  a  force  of  5600  N/m  (383  Ibf/ft) 
is  equivalent  to  1  1 .2  kPa  (1 .62  psi)  in  the  ice  cover. 
These  values  should  be  multiplied  by  a  factor  of  2.5 
to  take  into  account  the  variation  of  the  normalized 
load  in  the  experiments.  This  computation  results 
in  a  load  on  the  ICS  equal  to  14,000  N/m  (959  Ibf/ft) 
or  to  an  ice  pressure  ol  28  kPa  (4.06  psi)  lor  a  0.5-m- 
(20-in.-)  thick  icc  cover.  These  estimates  would  be 
valid  for  icc  cover  lengths  of  up  to  5  km. 

Ice  forces  on  ice  control  structure 
from  a  large  unconsolidated  ice  cover 

T'-c  s'  rcline  configuration  of  the  southern  end 
i  ‘  ^akc  Huron  is  sketched  in  Figure  13,  and  appears 
similar  to  that  considered  in  the  theoretical  model, 
r  rom  the  observations  made  in  the  field  and  the 
model  experiment,  we  can  expect  that  the  icc  cover 
the  ICS  will  be  trapped  between  the  shores  of 
■  'ike  and  the  ICS  in  the  form  of  wedge-shaped 
zones  (shown  as  shaded  areas  in  Fig.  13).  Probably 
these  zones  will  be  frozen  to  the  shore  bottom  due 
lo  formation  of  a  rubble  field  and  grounded  ridges, 
and  in  such  cases  any  force  transmitted  to  these 
zones  will  be  transmitted  to  the  shore.  On  the  other 
hand,  if  these  wedges  arc  not  frozen  but  remain  un¬ 
consolidated,  as  would  be  the  case  in  early  and  late 
winter,  the  ice  pressure  from  the  ice  cover  will  re¬ 
sult  in  pressure  on  the  ICS.  We  make  an  assumption 
that  the  icc  pressure  on  the  ICS  is  equal  to  the  icc 
pressure  exerted  on  the  stationary  ice  zones  from 
the  main  body  of  the  icc  cover,  i.e.  an  assumption 
of  an  isotropic  state  of  stress  in  the  wedge-shaped 
zone.  This  assumption  is  considered  to  be  conserva¬ 
tive. 

If  the  shear  stress  acting  on  the  ice  cover  is  large 
enough  to  create  a  state  of  stress  large  enough  for 
it  to  deform  in  its  plane,  a  slip  (or  shear)  line  will 
develop  between  the  stationary  icc  cover  adjacent 
to  the  ICS  and  the  main  body  of  the  ice  cover.  The 
boundaries  of  the  main  body  of  the  ice  cover  are 
approximately  straight,  and  we  can  make  use  of  the 
expression  developed  for  the  distribution  of  normal 
stress  components  along  the  boundary  (eq  A26). 

As  indicated  by  Figure  1 3,  we  have  assumed  that  the 
unconsolidated  ice  field  extends  to  about  50  km  and 
the  lower  end  of  the  icc  cover  is  at  1  km  from  the 
apex  of  the  converging  boundaries.  The  nondimen- 
sional  normal  stress  component,  ^00l^r2‘ri  ). in 
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Figure  13.  Sketch  of  southern  Lake  Huron. 


the  lower  end  of  the  ice  cover  is  approximately  equal 
to  0.08  (from  eq  A26  for  an  assumed  value  of  <t>  =  20° 
and  r2//-|  =  50).  Therefore,  Ngg  is  estimated  to  be 
5057  N/m  (347  Ibf/ft)  for  r2-rx  =  49  km  and  t  = 
1.26  N/m2.  Making  use  of  the  assumption  of  an 
“isotropic”  state  of  stress  in  the  wedge-shaped  zones, 
we  estimate  the  force  per  unit  length  of  the  ICS  to  be 
approximately  5000  N/m  (342  Ibf/ft),  which  is  equiv¬ 
alent  to  10  kPa  (1.45  psi)  of  pressure  in  a  0.5-m- 
(«=  20-in.-)  thick  ice  cover. 
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APPENDIX  A.  EQUATION  FOR  THE  STRESS  RESULTANTS  AND  VELOCITIES  OF  THE 
ICE  COVER 


The  purpose  of  this  appendix  is  to  present  the  equations  for  the  stress  resultants  and  velocities 
in  the  ice  cover  when  it  is  in  a  critical  state  (defined  by  the  state  of  stress  resultants  on  line  RS 
of  Fig.  2b)  and  to  derive  a  solution  of  these  equations  for  a  particular  case.  The  region  (r  =  r,  to 
r2  and  0  =  -6W  to  +0W)  of  the  ice  cover  is  shown  in  Figure  1 ,  and  the  ice  is  held  against  the 
boundaries  by  a  uniform  wind  stress  t  acting  on  the  ice  cover. 

The  governing  equations,  written  in  the  polar  coordinate  system  ( r ,  0)  shown  in  Figure  1,  are 
given  below.  In  these  equations  u  and  v  are  velocity  components  along  the  r  and  6  directions. 
Continuity: 


2"  +  "  +  1 
dr  r  r 


(Al) 


Equilibrium: 


,  1  3Ar0  |  A/rr  -  Ngg 
dr  +  r  dd  +  r 


=  T  cos  0 


dN'* 

dr 


1 


+  — 
r 


dNe_e:  2N^ 
dO  r 


=  -rsinfl  . 


(A2) 

(A3) 


Yield  criterion: 


+  Noo)  sin  -  Ngg) 2  +  4A&}1*  =  2c  cos0. 

Isotropy: 

o  _  /dv^  v  1  dtA/(~  9fA 
N,r  -Ng0  \dr  "  r  r  de)/\  drh 


(A4) 


(AS) 


where 

Nrr,  Nge,  Nrg  =  normal  and  shear  components  of  stress  resultants.  These  have  units  of  force  per 

unit  length  as  they  are  the  integrated  stress  across  the  thickness  of  the  ice  cover  and 
represent  the  force  interaction  between  floes  across  an  elemental  length  equal  to 
approximately  10  times  the  average  floe  size, 
r  =  the  uniform  wind  stress  per  unit  area  of  the  ice  cover  directed  toward  the  negative 
y  axis  as  shown  in  Figure  1. 

Ci  <t>  -  cohesive  strength  and  the  angle  of  internal  friction  of  the  ice  cover  during  in-plane 
shear  deformation  (see  Fig.  2  for  its  use  in  defining  the  yield  criterion). 


The  above  five  equations  are  to  be  solved  for  five  unknown  functions,  u,  v,  NJr,  Ngg  and  Nrg, 
which  must  also  satisfy  the  boundary  conditions.  The  solution  follows  a  procedure  similar  to  that 
of  Brennen  and  Pearce  (1978),  and  the  present  analysis  focuses  on  the  derivation  of  the  expressions 
for  the  stress  resultants  along  the  boundary  of  the  ice  cover  region.  Introducing  Sokolovski  functions 
(p,i£),  we  define  the  stress  resultants  in  the  following  way  such  that  the  yield  criterion  is  satisfied: 

A/rr  =  -p  +  (p  sin  <p  +  c  cos  <f>)  cos  2i^ 


^en  ~  -P  ~  ip  ^  c  cm  <j>)  cos  2 ^ 


Njo  =  [p  sin  <p  +  c  cos  6)  sin  2\fr 


(A6) 
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where  p  =  - (/Vrr  +  Nee)/ 2  and  i p  is  the  angle  between  the  greater  principal  stress  and  the  radial  di¬ 
rection. 

Using  eq  A6  and  introducing  a  =  p  +  c  cot  0,  we  may  rewrite  eq  A2,  A3,  and  A5  as  given  below: 


(1  -  sin  0  cos  20)  +  2o  sin  0  sin  2  C' 

or  or 

1  do  2  o 

-  sin  0  sin  20  -  rr  -  —  sin  0  cos  20  | 
r  oU  r 

(1  +  sin  ?  cos  20)  —  rr-  -  —  sin  0  sin  20 
r  OU  r 


-  sin  0  sin  2^  -2osincf  co: 
Or 

-3  u  3i'  v  1  3« 

2  r-  tan  2^  =  5-  -  -  +  -  zr  • 

3r  dr  r  r  30 


+  l)  +  T  COS  0  =  0 

*3  0  ) 

(A7) 

(°y-  + 1) 

\ur  / 

-  r  sin  0  =  0 

(A8) 

(A9) 

The  four  equations  (eq  A1,  A7,  A8  and  A9)  are  to  be  solved  for  four  unknown  functions,  u,  r,  o 
and  i^. 

It  is  to  be  observed  that  the  specific  values  ol  0  at  the  boundary  0  =  ±  9W  result  from  the  assump¬ 
tion  of  Coulomb  friction  on  the  boundary  of  the  form 

=  +  tan5  (A10) 

0  =  ±0W 


where  6  is  the  angle  of  friction  between  the  ice  cover  and  the  boundary.  From  the  eq  A6  this  leads 
to  the  relation 

IA"> 

Assuming  the  factor,  p!(p  +  c  cot  0),  is  approximately  equal  to  1 ,  we  may  evaluate  i^w  for  given 
values  of  8  and  0.  Further  assuming  that  the  velocity  and  stress  fields  are  symmetrical  with  respect 
to  the  centerline  0  =  0,  we  sec  that  0  increases  from  zero  at  U  =  0  to  the  known  value  of  0^  at  0  =  0W. 

Brennen  and  Pearce  (1978)  have  derived  a  solution  of  eq  A1,  A7,  A8  and  A9  by  expanding  u,  i; 
a  and  0  in  the  form  of  a  series  with  increasing  powers  of  9  We  shall  follow  the  same  procedure  to¬ 
gether  with  the  incorporation  of  the  simplifying  assumption  of  radial  flow,  i.e.  r  =  0,  and  linear  varia¬ 
tion  of  0  with  respect  to  0.  Thus  the  functions  u,  o  and  0  are  expanded  in  the  form 

W 

o=oQ(r)  +  o2(r)(jj2  +  oQ 


(A12) 

( A 1  3) 

(A14) 


The  procedure  is  to  expand  the  governing  equations  (eq  A1 ,  A7,  A8  and  A9)  in  powers  of  0  and 
to  obtain  ordinary  differential  equations  for  u0,  r/2,  o0  and  n-,  hy  equating  coefficients  of  powers 
of  0.  The  resulting  equations  arc  given  below: 


duo 

dr 


=  0 


(A15) 
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duQ 

u2  =  7^0v/r  — 


da0  2o0  /  \p  \ 

(1-sin0)  —  -  —  sin0(w^-)  +  r=° 


« *  P) 


(1  +  sin  0)  —  =  4  o0  0w  sin  0 


do0 

+  *K  sin  0  “^7"  +t/-0w  ■ 


The  solution  of  eq  A15  is  t/0  =  A/r,  and  the  function  u2  is  obtained  from  eq  A 16.  Using  eq  A12, 
the  radial  velocity  u  can  be  written  in  the  form 

(1  -20w  02/0w)  (A19) 

where  the  constant  A  is  indeterminate.  This  results  from  neglecting  the  inertial  terms  in  the  equa¬ 
tions  of  motion  (eq  A2  and  A3)  and  the  weak  coupling  between  the  stress  resultants  and  the  velocity 
field  through  eq  A5. 

The  solution  of  eq  A17  gives 

Trr'RT- 'Vo  («°> 

2  sin  0  (1  +  0W/0W) 

where  to  =  - -r. — : — -  and 

(I  -  sin  0) 

B  is  a  constant  to  be  determined  from  the  boundary  conditions.  The  expression  for  o2(r)  may  be 
obtained  by  substituting  a0(r)  in  eq  A1 8,  and  thus  deriving  an  expression  for  a(r,  (?)  from  eq  A1 3 
as  given  below: 

a  _  (fh )  T  (0\2  +  ^w)  sin  0  ^  02  (to  -  1 )  (1  -  sin  0)  |~j 

t/,  (to  -  1)  (1  -  sin  0)  I  \0W)  I  (1  +  sin  0)  2(1  +  sin  0)  IJ 


(r/r  i )  T  +/0\2  |0W.I36W  +  0W)  sin  < 

to  -  1 )  (1  -  sin  0)  \0W/  I  (1  +  sin  0) 

/rv-n  / 0\ 2 

Yl  /  [  +  \0j  1  (1  +  sin  0)  I 


A  different  expression  for  a  results  for  the  special  case  of  to  =  1. 

In  the  following,  the  expressions  for  the  stress  resultants  (/V^)  are  derived  for  the  special  case  of 
a  frictionless  boundary  condition  (i.e.  6=0  along  6  =  ±  0W  and  the  stress-free  boundary  condition 
(i.c.  /Vrr  =  0  and  Nt9  =  0)  along  r  =  r1  and  r2 . 

From  eq  A1 1 ,  we  get  0W  =  0  for  the  case  6=0,  and  thus  0  =  0  everywhere  in  this  region  as 
assumed  in  eq  A14.  When  0  =  0,  eq  A6  reduces  to  the  following: 

Nrr  -  c  cot  0  =  -(p  +  c  cot  0)  (1  -  sin  0) 

Ngg  -  c  cot  0  =  -(p  +  c  cot  0)  (1  +  sin  0) 


*r*=0. 


The  expression  for  a  =  (p  +  c  cot  <f>) is  derived  from  eq  A21  as  given  below: 


—  =  B  (—) 
Tri  Vi  / 


r\m  ^i  .  1  -  m 

r,  /  1-3  sin  0  2  1  +  m 
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(A23) 


2  sm  u 

where  /?;  =  - - -  . 

I  -  sm  0 


By  the  use  of  the  boundary  conditions  Nrr  (z, )  =  yVrr  (/-2 )  =  0,  the  constant  B  and  the  shear  stress 
rcan  be  evaluated  in  terms  of  the  geometrical  and  material  parameters.  The  shear  stress  t  necessary 
to  cause  a  critical  state  in  the  ice  cover  region  is  given  below: 


r(r2-r ,) 
C  cot  fS 


(1  -  m )/ 


(A24) 


The  stress  resultants  (/Vrr,  Ngg)  are  normalized  with  respect  to  r(r2  -z \ )  which  would  be  the 
maximum  pressure  in  the  ice  cover  region  if  its  behavior  were  similar  to  that  of  a  fluid  (i.e.  0  =  0, 
c  =  0)  and  the  ice  cover  were  blocked  at  z  =  z,  (i.e.  u  =  0  at  z  =  z, ).  The  normalized  expressions  for 
the  stress  resultant  are  given  below: 


No<>  /  1  \  j  1_-_m  dfl  prw  +  1)z-z,  (m  +  1 )  rm -rfl 

r(r 2~r])  \1  -m)  1  +  m  2J  r2  -z,  J 


(A26) 


The  graphs  of  Nlr/T(r2  -r, )  and  N60/T(r2  -z,)  are  given  in  Figure  Ala  with  respect  to  z/z,  for 
various  values  of  0/6w  and  for  given  values  of  the  angle  of  internal  friction  6  and  r2/r]  ratio.  The 
distiibution  of  stress  resultants  along  the  boundary  0  =  t  0W  has  been  plotted  in  Figure  A1  b  for 
various  values  of  the  angle  of  internal  friction  6  and  for  a  given  value  of  thez2/z,  ratio. 

A  simplified  analysis,  similar  to  that  by  Reimer  et  al.  (1979),  is  given  below.  In  this  analysis  the 
following  assumptions  have  been  made:  1 )  the  unknown  variables  are  independent  of  0,  2)  the  flow 
is  radial  (i.e.  v  =  0),  3)  the  stress  resultants  Nu  and  Ngg  are  the  principal  stresses  (i.e.  Nrf)  =  0),  and 
4)  the  shear  stress  facts  radially  inward.  Equations  A1-A5  then  reduce  to  the  following  set  of  equa¬ 
tions: 


Continuity: 


du  u_ 
dr  r 


=  0 


Equilibrium: 

dNr,  Nlr  -  Nu0 

—  + - =  r  . 

dr  r 

Yield  criterion: 

Nn  -  c  cot  0  _  sjn  0 
Nnf)  -  c  cot  0  1  +  sin  0 


(A27) 


(A28) 


(A29) 


The  solution  of  eq  A27  is  u  =  A/r  where  the  constant  of  integration  A  is  indeterminate. 
Eliminating  Nnu  from  eq  A28  with  the  help  of  eq  A29  we  get 


dr 


cot  4>)  -  m  (Nn  -c  cot  0)/r  =  r 


(A30) 


where  m  =  2  sin  0/(1  -  sin  0).  The  solution  of  eq  A30  may  be  written  in  the  following  form: 


Nn  -  c  cot  0 
Trt 


.  nr , 

(1  -m) 


(A31) 


By  the  use  of  the  boundary  conditions  /Vfr(r, )  =  Nrr(r2)  =  0,  we  may  evaluate  B  and  the  shear  stress 
t  in  terms  of  the  material  and  geometrical  parameters  of  the  ice  cover  region.  The  shear  stress  neces¬ 
sary  to  cause  a  critical  state  in  the  ice  region  is  given  by 


f{r 2  ~ri) 
c  cot  0 


=  (1  -m)l 


(A32) 


The  above  expression  is  identical  to  those  given  by  Sodhi  (1975)  and  Reimer  et  al.  (1979).  The 
expressions  for  Nrt  and  N00  arc  normalized  with  respect  to  r(r2  -r, )  and  are  given  below: 


Nr 


T{r2-r. 


N. 


oo 


f(r t  -r , 


(T 


1  [r-r i  r™  -ryl 

-m)  [r2-r,  rm  .  rm  J 

1  p1  +  m)  r  -  r,  (i  +  m)  r">  -  r™1 

■rt?)l  r2-r,  rm  „  rm  I 

u  r2  '  1  J 


(A33) 


(A34) 


The  above  expressions  are  identical  to  eq  A25  and  A26  except  for  a  factor  involving  an  expression 
of  0. 
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APPENDIX  B.  MONTHLY  SUMMARY  OF  WIND  DATA  AT  PORT  HURON. 


WSW  7.1  50  SO  0  0  0  0 


February  1980 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Sodhi,  D.S. 
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available  from  National  Technical  Information  Ser¬ 
vice,  1982. 
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